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Abstract 

Self-compression of multi-millijoule femtosecond laser pulses and dramatic increase of the peak intensity 
are found in pressurized helium and neon within a range of intensity in which the ionization modification of 
the material parameters by the pulse is negligible. The pulse propagation is studied by the (3+l)-dimensional 
nonlinear Schrodinger equation including basic lowest order optical processes - diffraction, group velocity 
dispersion of second order, and Kerr nonlinearity of third order Smooth and well controllable pulse propa- 
gation dynamics is found. Constructing of compressed pulse of controllable parameters at given space target 
point can be achieved by a proper chose of the pulse energy and/or gas pressure. 

1 Introduction 

The propagation of high-intensity femtosecond laser pulses in a bulk medium is a complex phenomenon accom- 
panied by a dramatic spatiotemporal and spectral reaiTangement of the pulse. This originates from the highly 
intensive nonlinear and dispersive interaction of the optical pulse with the medium. The rigorous description 
of the pulse propagation is based on the (3h-1)D nonlinear Schrodinger equation (NLSE), in which diffraction, 
group velocity dispersion (GVD) of second order, and Kerr nonlinearity of third order form a basic set of optical 
processes of lowest order [1-6], which is universal for all media and parameters of the input pulse and belongs 
to the strongest effects acting on the pulse in a centrosymmetric neutral medium. Additional terms in the NLSE 
are required to describe some particular features or more extreme conditions in the pulse propagation as, e.g., 
self-steepening, space-time focusing, non-instantaneous nonlinearity, higher order dispersion, etc. [7-9]. Higher 
orders of nonlinearity and ionization must also be considered for pulses creating significant ionization [10-12]. 
While the femtosecond pulse behavior strongly depends on the material and pulse parameters, some general 
simple understandings were widely accepted. In this way, simultaneous collapse of the pulse in space and time 
and formation of a spatiotemporal soliton is predicted at negative GVD [1 1, while, in the more common case 
of positive GVD, the theoretical and experimental studies predict time broadening and splitting of the pulse 
[2-5,7,8,10]. In a single study, initial pulse shortening is found while not supported by experimental evidences, 
instead, splitting of the femtosecond pulse is shown |j6J . 

Self-compression (SC), before splitting, of high-intensity femtosecond laser pulses in positive GVD has 
been discovered in various types of media - pure atomic and molecular gases [13-15] and fused silica lfT6l . 
This allows recognizing the SC as a new phenomenon - a kind of a time focusing, or the temporal counterpart 
of the usual space domain self-focusing (SF) lfT4l . The mechanism of both phenomena is different but, as 
will be shown below, the SC is closely related to the SF. The generation of compressed (almost five times in 
magnitude) pulses in the experiments in gaseous media is accompanied by a strong increase of the peak intensity, 
an improvement of the spatiotemporal pulse shape, and a stable propagation of the compressed pulse along a 
distance that is many times longer than the characteristic length of the strongest pulse rearrangement factor, 
in this case - the nonlinearity [14j . The time-spectral studies show that the compressed pulses are transform- 
limited ifTTll . SC of high-intensity femtosecond laser pulses was also observed in air ifTSl , and in ionized noble 
gasses [19-22]. Recently, SC was also reported in BK-7 glass using negatively chirped pulses f23). Generation 
of long light filament and twofold pulse compression has been observed experimentally in air [181 . Substantially 
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stronger pulse compression has been achieved in ionized noble gasses [19-22]. Generation of few-cycle pulses 
through filamentation is achieved experimentally in argon- filled gas cell ||20l . Generation of single-cycle pulse 
is predicted numerically in argon gas with pressure gradient f2Tl. Grate attention has been paid in the study of 
the femtosecond pulse splitting. From the other side, while the SC is much more interesting for practical use, 
its exact "physical mechanism ... is still not fully understood" Il22ll . Recently, we have shown for the first time 
that SC and intensity gain of high intensity femtosecond laser pulses, before the pulse splitting, can be achieved 
in a positive GVD medium based on basic set of the lowest order optical processes, diffraction, group velocity 
dispersion of second order and Kerr nonlinearity of third order, providing the GVD of the medium is low enough 
||24|. The physical mechanism of SC in this case was found and referred to as SC in low dispersion regime. 
The numerical simulations were performed for pressurized argon assuming sub-millijoule femtosecond laser 
pulses so as to avoid the effect of the ionization on the pulse propagation. The lack of ionization substantially 
simplifies the pulse propagation dynamics. To scale up the above method to higher energy/intensity pulses, 
media of higher ionization potential have to be used. Here we demonstrate SC of multi-millijoule femtosecond 
laser pulses in pressurized helium and neon. The pulse compression is accompanied by a dramatic increase 
of the peak intensity at almost no loss of energy. The pressurized helium and neon possess the lowest specific 
GVD among the known media, which provides an increased degree of SC and peak intensity gain in comparison 
with argon. The dynamics of the optical pulse at lack of ionization is relatively smooth and it allows a reliable 
control on the pulse propagation and constructing a pulse of well predictable parameters in given space target 
point. The degree of SC and intensity gain, as well as the space positions of maximal SC and intensity gain of 
the pulse can be controlled by the pulse energy and/or gas number density/pressure. 



2 Physical Model 

We consider a physical model of the high-intensity femtosecond pulse propagation in a centrosymmetric non- 
Unear medium, which includes the basic set of lowest order optical processes, diffraction, GVD of second order 
and Kerr nonlinearity of third order The NLSE (in moving frame) 
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for the complex field amplitude E{r, z, r) of a pulse propagating along the z-direction, including the above 
mentioned physical processes, is used as a propagation equation within that model. Here, k — 27r/A is the 
magnitude of the wave-vector, n is the linear refractive index of the medium at the central wavelength A of the 
pulse, is the transversal Laplacian, (32 is the GVD, and 712 is the nonlinear refractive index. The field is 
normalized such that / ~ \E\'^ is the intensity. Neon and helium are used as nonlinear media. The nonlinear 
refractive indexes of neon and helium, 712 (Ne) = 7.8 x 10^^^ Pcm-^/W and n2(He) = 4.2 x 10^^^ Pcm^AV, 
have been determined from a comparative study of the relative nonlinear refractive indexes of noble gas media 
II25J . taking the nonlinear refractive index of Ar, 712 (Ar) = 9.8 x 10^^° Pcm^AV, as a reference value. The 
GVD of neon and helium has been determined from the Sellmeier equation ll26l and found to be P2 (Ne) =0.02 
Pfs^/cm and /32(He) = 0.01 Pfs^/cm. In the above expressions P is the pressure in atm. The initial pulse is a 
linearly polarized chirp-free Gaussian E{r, z = 0,t) = Eq exp(— r^/2rg — t'^ /2tq) in space and time having 
axial symmetry and 100 fs time duration (full width at half maximum (FWHM)) of the intensity profile. The 
central wavelength of the pulse is tuned at 800 nm. 

The proper initial conditions are important in order to remain within the range of validity of the present 
model. For such a purpose, special attention is paid to the ionization. The ionization strongly and non- 
instantaneously modifies the material parameters and thus, the pulse propagation. Also, it is highly nonlinear 
process, which means that it is very sensitive to the changes of the intensity/field strength and its control is 
very difficult in practical means, more over - the pulse intensity strongly changes during the propagation. The 
avoiding the ionization of medium can be achieved by a suitable choice of the input pulse energy and/or the 
gas number density/pressure so that the peak intensity, raised by the SF, does not surpass the level at which the 
ionization causes substantial modification of the material parameters. For instance, the plasma contribution to 
the refractive index should be kept negligible in comparison to the respective contribution of the cubic nonlinear 
Kerr effect within the whole range of intensity variation. Avoiding the ionization ensures a relatively simple, 
smooth and, thus, well controllable pulse propagation dynamics. 

The determination of the absolute value of the ionization rate and, in this way, the plasma density created 
by the pulse itself is difficult task. The description of the ionization of atoms and simple molecules by the 
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general ionization theories is not quite satisfactory [27J . The ionization rate and the plasma number density pre- 
dicted, e.g., by Perelomov-Popov-Terent'ev (PPT) theory ||28l . which appears to be among the most successful 
ionization theories [27 1, closely matches the trend of the experimental data over a wide range of intensities - 
from multi-photon to tunnel ionization regimes. However, the agreement between the experimental data and the 
theoretically predicted plasma number density, based on the PPT theory, is still within a factor of two. The dis- 
crepancy with the other ionization theories may reach orders of magnitude ll27l . From the other side, precision 
measurements of the ionization of helium [ 29 J as well as of the other noble gases ||27l are performed over a wide 
counting range. Unfortunately, only the relative value of the ion signal has been linked to the absolute value of 
the peak intensity of the pulses in these experiments [27,29], whereas, due to the ion detection technique, the 
absolute value of the plasma density remains undetermined. To account for the latter, we require at least one 
absolute reference point to calibrate the relative plasma number density-peak intensity relationships found from 
these experimental studies. Then, the precise relative dependence can be used to determine the respective abso- 
lute value dependence. To find such a reference point, the following approach will be applied. The intensity of 
complete ionization of neon and helium atoms will be determined using the simple barrier suppression ioniza- 
tion (BSI) model f3U\. Within the BSI model, the threshold intensity /bsi at which the optical field suppresses 
the Coulomb potential in an atom so as to allow a bound electron to escape freely above the potential barrier 
(above barrier ionization) without tunneling is (in atomic units) 

-^BSI = Y^' (2) 

where Ip is the ionization potential of the initial atom/ion, and Z is the charge of the created ion. At the ion- 
ization potential of neon, Ip = 21.56 eV, and helium, Ip = 24.58 eV, the intensity threshold for complete 
ionization of the medium is found to be /Bsi(Ne) = 0.86 x 10^^ W/cm^ and /Bsi(He) = 1.46 x 10^^ W/cm^, 
respectively. The BSI intensity of helium atom so determined is relatively close to the experimentally deter- 
mined value of 0.8 x 10^^ W/cm^, beyond which indication of saturation of the ionization can be distinguished 
(the slow /■^/^ increase of the experimentally measured ion yield beyond the sited intensity is attributed to the 
expansion of the Gaussian focal volume) ||29l . The intensity of the barrier suppression ionization, /bsi. cor- 
responds to complete ionization of the medium, and, at / = /bsi. the plasma number density will be equal 
to the number density of the initial neutral gas (the plasma recombination will be neglected on the time scale 
comparable to the pulse duration of 100 fs). The latter is, in fact, the reference point that will be used here in 
order to find the relationship between the absolute values of the plasma number density and the peak intensity 
of the pulse. Using such a relationship allows to account for the intensity at which the plasma contribution to 
the refractive index is still negligible. 

According to the above considerations, the electron/plasma number density A^e at / = /bsi will be equal to 
the initial number density of the neutral gas, Nno- We will take tentatively 6 orders of magnitude lower plasma 
number density than Nno, and the intensity at which such plasma number density is created will be determined 
approximately from the respective experimental curves. Figure 1 of Ref.[29[ and Figure 5 of Ref.[27[ for He, 
and Figure 4 of Ref.[27[ for Ne. From the sited figures, the intensity corresponding to such a plasma density 
was found to be about I = 2 x 10^^ W/cm^ and / = 1 x 10^* W/cm^ for helium and neon, respectively. 
At the intensity so obtained and the respective plasma number density, the magnitude of, otherwise, negative 
plasma contribution to the refractive index, Aupi « ujpi/2uj'^ = ATie^Ne/lmu}'^ [31], will be compared to the 
respective Kerr effect contribution to the refractive index, n„; ~ n2\E\'^. This is justified taking into account 
that the latter quantities determine the magnitude of the contributions of these effects to the respective terms, 
k{ujli/uj'^)E and kn2\E\^E, in the NLSE [12[. For 30 atm gas pressure of helium (the maximal gas pressure 
used here), the plasma contribution to the refractive index at / = 2 x 10^"* W/cm^ {i.e., assuming, as sited 
above, 6 orders of magnitude lower plasma density than Nno) is ripi « us^ij^uP' = 2.15 x 10~^ while the 
Kerr effect contribution to the refractive index at the same intensity is Uni = n2\E\'^ = 2.52 x 10~^. The same 
quantities for 17 atm pressure (around which most of the simulations are performed here) are n^i = 1.22 x 10^^ 
and Uni = 1.43 x 10^^ for the case of heHum at / = 2 x 10^** W/cm^ intensity, and Upi = 1.22 x 10^^ and 
Uni = 1.33 X 10~^ for the case of neon at / = 1 x lO^'' W/cm^ intensity. As can be seen, the plasma contribution 
to the modification of the refractive index is more than 2 orders of magnitude lower than the respective Ken- 
effect contribution and neglecting the plasma contribution seems justified for neon and helium at intensities not 
exceeding / = 1 x lO^"' W/cm^ and I — 2 x 10^^ W/cm^, respectively. 

Within the above mentioned intensity range we are still in perturbative regime what concerns the neon and 
the helium atoms. The relation between two successive terms of the nonlinear polarization for bound-bound 
transitions is given by the following approximate relation [32J abb = ~ {eEaab / hAuj)^ 
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(for an isotropic medium the even order nonlinearities are zero), where Ea is the time dependent field amplitude, 
flb is the Bohr radius, and TiAoj = hueg — hu) is the energy detuning of the laser photon energy Tiu) (at the carrier 
frequency of the pulse) from the atomic resonance of energy huoeg between the nearest exited atomic state 
(coupled to the ground state by a nonzero dipole moment) and the ground state. For the helium atom, the lowest 
excited level with nonzero dipole moment to the singlet ground state Is is the singlet 2p state, lying about 21.2 
eV above the ground state. At the field strength corresponding to the maximal intensity of / = 2 x 10^'' W/cm^, 
the relation between the leading cubic term and the next quintic term is about af,b ~ 10^^, or the truncation of 
the nonlinearity after the cubic term seems justified within the present considerations. One approach to control 
the peak intensity below the sited levels is to change the gas pressure while keeping the input energy fixed for 
given study. This is also a convenient approach for a real experiment. 

In view of the above considerations, the propagation equation ([1) (which does not include irreversible losses 
and non-instantaneous processes) is appropriate to describe the propagation of femtosecond pulses within the 
conditions in which the present model is valid. The lack of substantial losses allows imposing energy conserva- 
tion condition ll24l 



where Wq is the input pulse energy, W{z) is the pulse energy along the propagation distance, and the integration 
is taken over the transversal cross-section s and the local time r of the pulse. The latter also represents a 
normalization condition for the complex field envelope E along the propagation distance z. Applying the 
energy normahzation condition for the field envelope allows to perform absolute value calculations. 

3 Results and Discussion 

The NLSE ^ is solved numerically in a polar coordinate system (assuming axial symmetry) by using a split- 
step method similar to that applied in ll24l . Nonlinear media of highest ionization potentials - pressurized neon 
and helium, are considered in this study. For better comparison of the pulse compression capabilities of neon 
and helium, completely identical initial conditions for the optical pulse and the medium have been taken for the 
initial numerical simulations in both media. Our simulations show that at Wq = 2 mJ input pulse energy and 
P = 17 atm pressure of neon or helium, the maximal peak intensity does not surpass / = 1 x 10^* W/cm^ in 
neon and / = 2 x 10^"* W/cm^ in helium along the whole propagation distance considered here. Consequently, 
at such initial conditions, the parameters of the pulse fall well within the range of validity of the present model. 

The SC in low dispersion regime is inherently accompanied by a gain of the peak intensity of the pulse. 
That is why, two parameters have main importance here - the magnitude of the peak intensity gain fcic and the 
magnitude of the SC, ksc- The peak intensity gain will be defined as the relation between the peak intensity 
of the pulse at given space point z, Io{z), and the peak intensity of the input pulse, /o(2: — 0), i.e., kiQ{z) = 
Io{z)/Iq{z = 0). The magnitude of the SC will be defined as the relation between the time duration of the 
pulse (FWHM) at given space point z, tq{z), and the time duration of the input pulse to(2: = 0), i.e., ksc{z) = 
Toiz = 0)/to{z). 

The evolution of the transversal pulse width 2ro (FWHM at r = plane), the normalized peak intensity 
\Eo{r = 0, z,T — 0)/ Eo{r = 0, z = 0,T = 0) (i.e., the peak intensity gain), and the pulse duration tq (FWHM) 
of the intensity profile along the propagation distance z at 2 mJ input pulse energy and 17 atm pressure are 
shown in Figures l(a),(b),(c) for neon and helium. The pulse shows similar behavior in both media, which is 
also close to the qualitative behavior of the pulse in the case of argon [24|. The pulse rearrangement begins 
with a confinement of the pulse in the transversal direction due to the SF, Figure 1(a). The minimal width of 
the pulse in the transversal direction is reached at z = zgp = 0.653 m in neon and z = zgF = 1101 m in 
helium. The SF confines the peripheral energy of the pulse toward the longitudinal axis, which leads to an 
increase of the peak intensity. Figure 1(b). The space positions at which the peak intensity reaches its maximal 
value is ziq = 0.642 m in neon and zjq — 1.089 m in helium. The corresponding maximal peak intensity 
gain at these positions in neon and helium is approximately kiQ{z — ziq) = 29 and kiQ{z — ziq) — 48 times, 
respectively. The most interesting and nonintuitive result is that the pulse undergoes SC in time (before splitting) 
while propagating along a positive GVD medium. Figure 1(c). The position of the maximal SC is zsc = 0.620 
m in neon and zsc = 1-067 m in helium. The magnitude of the SC in these positions is fcsc(^ = zsc) = 3.4 
times in neon and ksc{z = ^sc) — 4.4 times in helium. The intensity gain at the position of maximal SC is 
approximately kiG{z — zsc) — 24 times in neon and kiQ{z — z^c) = 37 times in helium. As can be seen, the 
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positions of the maximal SC and intensity gain do not coincide but are relatively close. This means that within 
given propagation distance simultaneous time compression and intensity gain of the pulse can be achieved. In 
a "chronological" order, first appears the maximal SC at zsc, next occurs the maximal intensity gain at ziq, 
and finally - the maximal self-focusing at zsf- Once the pulse reaches the maximal SC at z = zsc, it begins 
broadening in time from zgc to ziq and the magnitude of the SC reduces to ksc{z — ziq) — 2.8 times in 
neon and fcscl^ — ^ig) = 4 times in helium. At the same time, the intensity continues to rise until z — ziq 
because the pulse continue to narrow in the transversal direction due to the SF while it reaches the position of 
maximal SF at z = zsf- The propagation of the pulse beyond ziq leads to the usual pulse splitting [2-8,10,24]. 
Indication of the pulse splitting in Figure 1(b) is the rapid fall of the peak intensity at the central peak point of 
the pulse at r = beyond z = zig- 
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Figure 1 : Evolution of the transversal pulse width (FWHM) (a), normalized peak intensity (b), and pulse du- 
ration (FWHM) (c) versus the propagation distance. The positions where the respective quantities are not well 
defined because of the pulse splitting are indicated with dashed lines in the curves. 



The comparative studies of the pulse compression in neon and helium clearly reveals a well expressed 
trend — the pulse propagation dynamics develops more slowly in helium than in neon, while the magnitude 
of the pulse compression and the intensity gain is substantially higher in helium than in neon at equal other 
conditions. This can be explained by the lower GVD of helium. The latter diminishes the effect of dispersion 
spreading of the pulse while it effectively compresses in time due to the SF. In fact, helium possesses the lowest 
specific GVD (at equal other conditions) and the highest ionization potential among the known media. That is 
why it appears to be the top "speed" medium within the present concept of propagation type compression of 
high intensity femtosecond laser pulses. 

The optical pulse do not tend to form light filaments within the basic set of processes, but starts to diverge 
once it reaches the position of maximal SF, z — zsf. Figure 1(a). That is why, one may speculate that the 
generation of light filaments results from the defocusing effect of the above cubic nonlinearity, e.g., x''^\ ™d 
the ionization. These processes provide an instantaneous dynamical (due to the x*^^^ -nonlinearity) and non- 
instantaneous (due to the ionization) balance of the x'^'^-'-SF and lead to formation of more or less stable light 
filaments [14,18,20]. For the experimental conditions considered here, the intensity raised by the x^'^^-SF does 
not reach high enough values so as to "switch on" the x*^^-* and the ionization defocusing, and the x'^^-SF is 
restricted only by the diffraction and the pulse splitting due to the dispersion. The latter factors seem not to be 
sufficient to form relatively stable light filaments. 

A general explanation of the SC in low dispersion regime has been found ll24l based on the energy conserva- 
tion. For the case of localized electromagnetic fields (pulses), the energy conservation has a simple geometrical 
sense, expressed by Eq.(l3]i. According to it, the energy integral can be considered as a four-dimensional "vol- 
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ume" (two transversal directions over the cross-section of the pulse, a longitudinal direction - the local time, and 
a "vertical" direction - the intensity \E\'^), which must be conserved along the propagation coordinate z. An in- 
tensive rearrangement of the pulse is triggered by the SF coming from the nonlinear term. The rapid (relatively 
to the dispersive broadening) development of the SF is the main point in the pulse compression mechanism 
considered here. The SF confines the pulse radially from all transversal directions toward the longitudinal axis. 
This increases the peak intensity and further accelerates the pulse rearrangement. The increased peak intensity 
intensifies the self-phase modulation thus generating new spectral components - the "red"-frequency compo- 
nents are generated on the raising edge and the "blue"-frequency components are generated on the falling edge 
of the pulse. In a positive GVD medium, this results in faster propagation of the fore edge and slower propa- 
gation of back edge of the pulse, which tends to cause pulse broadening. Due to the low dispersion of the rare 
gases neon and helium, the pulse would not substantially expand in time (within given propagation distance) 
if it is ruled solely on the dispersion. At the same time, the strong confinement of the pulse in the transversal 
direction due to the SF results in a strong increase of the pulse intensity. Figure 1(b). The latter leads to an 
effective, but real, shortening of the pulse (based on the standard FWHM characterization of the pulse duration) 
before the pulse succeed to broaden substantially due to the dispersion. This is the mechanism of the pulse com- 
pression observed here. Such a concept of femtosecond pulse compression, formulated for the first time in ll24l . 
has been called self-compression in low dispersion regime. The pulse compression in low dispersion regime is 
inherently accompanied by a strong increase of the peak intensity at almost no loss of energy, in accordance to 
the condition of conservation of the pulse energy. This is the main advantage of the present pulse compression 
method. In fact, the pulse behavior is subject of competition between two processes acting in opposite direc- 
tions - the pulse compression due to the intensity gain as a result of the SF and the pulse broadening due to 
the dispersion. The final outcome depends on the relative rates of both trends. Pulse compression along given 
propagation distance can be achieved if the rate of pulse broadening due to dispersion remains lower that the 
rate of pulse shortening due to intensity gain forced by the SF. Otherwise, pulse broadening and splitting will 
be observed. This has been confirmed numerically |24|, showing that if the GVD exceeds given value (at equal 
other conditions), only pulse broadening and splitting, but not SC, is observed. This means that the SF does not 
necessarily lead to pulse compression [33,34]. The pulse compression in low dispersion regime is substantially 
related with the low value of GVD and the efficiency of this mechanism increases when the GVD diminishes. 
This explains why the degree of SC and the intensity gain increases from argon ll24l toward neon and helium. 
Figure 1. It is important to underline that the above SC phenomenon is substantially (3h-1)D effect that is based 
on the SF and cannot be observed in the (1h-1)D case at positive GVD {e.g., optical fibers), where the SF is 
not considered. The latter also shows the mutual relation between the different dimensions in the high-intensity 
femtosecond pulse dynamics. 

Although the pulse collapses in all dimensions in space and time, the confinement mechanism of the pulse 
within the basic set of optical processes is still incomplete. This is because confinement "force" in this case, 
i.e., the SF, exists only in the transversal space dimensions but not in the time dimension at positive GVD. The 
SC in the time dimension results from the retarded dispersion broadening of the pulse due to the low GVD of 
the medium. Consequently, the stable pulse propagation in the positive GVD medium, as has been observed 
experimentally fT4], requires additional process(es) to stabilize the optical pulse. 

The present mechanism of pulse compression is based on the geometrical rearrangement of the pulse, in 
which the nonlinearity by means of the SF plays a leading active role, whereas the GVD, while it is also very 
important, has a secondary role, i.e., to promote or to hamper the SC. That is why, such kind of SC is expected 
to work not only for initially chirp-free pulses, but for (regularly) chirped pulses and pulses having complicated, 
even irregular, frequency distortions. In the latter cases, the efficiency of such SC mechanism is subject of 
additional studies. The specified basic set of optical processes is universal for all centrosymmetric media and 
experimental conditions. In particular, the basic lowest order processes are common for cases without or with 
ionization. That is why, the mechanism of SC described here should also play important role for the cases with 
substantial ionization of the medium. To what extend the proposed here mechanism of SC plays important role 
for the cases beyond the basic set, besides on the particular experimental conditions, depends also on existing of 
other, if any, accompanying process(es) of comparatively strong action as those from the basic set. Among such 
additional processes, the most important are expected to be the ionization and above-cubic nonlinearities fT?). 
The contribution of the ionization and the above-cubic nonlinearities (in the case of higher intensity pulses) 
to the present mechanism, as well as the efficiency of the latter, is subject of additional studies. No principle 
limitations seem to exist for the present mechanism of pulse compression to work for much shorter pulses then 
the ones used here but more advanced propagation equation [32 1 have to be used in that case. 

A kind of time shortening can also be found within the pulse splitting phenomenon, even if the initial pulse 
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has newer reached shortening but broadens continuously and finally - splits. At given stage of development 
of the split pulse, each of two sub-pulses, in which the initial pulse splits, may become shorter than the initial 
pulse. The split pulse, however, is usually considered as deteriorated pulse, having reduced energy and peak 
intensity in each subpulse, and such a kind of compression seems useless for the practical application. The SC 
before splitting, as it is observed here, must be clearly distinguished from the derivation of shorter subpulses 
within the pulse splitting phenomenon. 
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Figure 2: Dependence of the maximal peak intensity and the respective peak intensity gain versus the input 
energy of the pulse at fixed pressure of helium. 
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Figure 3: Dependence of the space position ziq, where pulses of maximal intensity are generated, versus the 
input energy of the pulse at fixed 17 atm pressure of helium. 



For comparison with the results for 2 mJ energy pulses, compression of higher energy 5 mJ-pulses were also 
studied in neon and helium, keeping the other conditions same as above. In this case, we found lower degree 
of the peak intensity gain, kiQ{z = ziq) = 13 times in neon and kiQ{z = ziq) = 22 times in helium, which 
results in lower degree of pulse compression, ksc{z = zsc) = 2.3 times in neon and fcsc(^ = -^sc) = 3.0 
times in helium. Whereas the initial peak intensity is higher for 5 mJ input energy, the maximal peak intensity 
at ZiG does not substantially increases when the pulse energy increases from 2 mJ to 5 mJ. Thus, the maximal 
peak intensity at ziq for the case of 2 mJ input pulse energy is 1 x 10^^ W/cm^ in neon and 1.9 x lO^"* W/cm^ 
in helium, while for the case of 5mJ these values are 1.2 x 10^^ W/cm^ and 2 x 10^'' W/cm^, respectively. 
The pulse compression is also lower at higher input energy. Consequently, the simple increase of the pulse 
energy does not necessarily lead to increase of the SC and the intensity gain and a proper optimization of the 
parameters of the initial pulse and the medium is required in order to achieve the highest performance from the 
pulse rearrangement process. In view of the non-intuitive SC behavior of the high-intensity femtosecond pulses 
versus the input pulse energy and medium density, a more detailed study on that subject is required and will be 
presented below. 
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Whereas the propagation of high intensity femtosecond laser pulses in the heavier rare gases (argon, krypton 
and xenon) has been widely investigated [10,11, 13-15, 19-22,24], the lighter rare gases (neon and, especially - 
helium) remain insufficiently explored as propagation media. From the other side, the above results show the 
superior self-compression and intensity gain capabilities of helium in comparison with that of neon and argon 
Il24l within the present pulse compression concept. That is why, the further studies will be focused namely on 
helium. In the following, the dependence of the main parameters of the pulse on the input pulse energy and 
helium pressure will be shown. The development of the maximal peak intensity at z = ziq and the respective 
magnitude of the peak intensity gain versus the input pulse energy at fixed helium pressure, P = 17 atm, are 
shown in Figure 2. As one can see, the absolute value of the peak intensity increases with the input pulse 
energy, while, at the same time, the magnitude of the peak intensity gain, kiQ, decreases. The increase of the 
peak intensity with the pulse energy is more rapid in the low energy scale and tends to saturate at high energies. 
The space position, zig, where a pulse of the maximal peak intensity is created, as a function of the same values 
of the input pulse energies, as in Figure 2, is shown in Figure 3. The space position ziq, where the pulse has 
maximal intensity gain, changes in a kind of reciprocal way with the pulse energy. At relatively low energy, 
say Wq = 1 mJ, the position of the maximal peak intensity is very far from the entrance of the medium, and it 
rapidly approaches the entrance when the pulse energy increases. Figure 3. 
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Figure 4: Dependence of the minimal pulse duration and the respective coefficient of SC versus the input energy 
of the pulse at fixed 17 atm pressure of helium. 
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Figure 5: Dependence of the space position zgc, where pulses of minimal pulse duration are generated, versus 
the input energy of the pulse at fixed 17 atm pressure of helium. 



The dependence of the minimal pulse duration at z = zgc and the respective magnitude of the pulse com- 
pression versus the input pulse energy at fixed helium pressure, P — 17 atm, are shown in Figure 4. The space 
positions, zsc, where a pulse of minimal pulse duration is created, at the same input pulse energies, as in Figure 
4,is shown in Figure 5. In contrast to the pulse intensity. Figure 2, the minimal pulse duration increases almost 
linearly with the input pulse energy. What is also interesting, the degree of the SC, ksc, decreases with the in- 
crease of the pulse energy. At the same time, the positions zsc where maximally compressed pulses are created. 
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change in a similar reciprocal way in Figure 5 as for the case of ziq in Figure 3. At the lowest energy used, 
Wq = 1 mJ, zsc is very far from the entrance of the medium and it very rapidly moves toward the entrance 
when the energy increases. Figure 5. 
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Figure 6: Dependence of the maximal peak intensity and the position ziq, where such pulses are generated, 
versus the gas pressure of the medium at fixed Wq = 5 mJ input energy of the pulse. The coefficients of the 
intensity gain fcic at the respective values of gas pressure are also given in the figure. 
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Figure 7: Dependence of the minimal pulse duration and the position zjg, where such pulses are generated, 
versus the gas pressure of the medium at fixed Wo = 5 mJ input energy of the pulse. The coefficients of the 
self-compression ksc at respective values of gas pressure are also given in the figure. 



The dependence of the maximal peak intensity al z — ziq as well as the position ziq, where such pulses 
are generated, from the helium pressure at fixed input pulse energy of Wq = 5 mJ is shown in Figure 6. The 
dependence of the minimal pulse duration at z = zsc as well as the position z^c, where such pulses are 
generated, on the helium pressure at fixed input pulse energy of Wq = 5 mJ is shown Figure 7. The coefficients 
of the intensity gain, fcjG, and the SC, fcsc, at the respective gas pressure, are also given in the figures. These 
studies show that the strongest degree of SC and intensity gain is achieved at lowest pressure, P = 15 atm, 
which also means lowest value of the total GVD of the medium. The peak intensity at the lowest helium 
pressure in Figure 6, as well as at the highest energy pulses in Figure 3, surpasses the accepted here value of 
maximal intensity of / = 2 x 10^'' W/cm^ (for the case of helium) within about 10%, and these results are shown 
for completeness of the study. Further reduction of the pressure results in peak intensity which is substantially 
beyond the specified value, and, for Wq = 5 mJ pulses, the simulations were not performed below that pressure. 

The evolution of the spatiotemporal pulse shape with the variations of gas pressure at fixed Wq = 5 mJ 
input pulse energy has been also investigated. A "time view" of the spatiotemporal pulse shape taken at the 
position of maximal SC, z — zsc, is shown in Figure 8. The parameters of the predicted pulses at the specified 
space points are also given in the figure. These results demonstrate the rearrangement capabilities of the present 
pulse formation method at strong nonlinear interaction with the medium and low dispersion regime. If the 
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medium ends around z^c ov within zsc ^ ^IG interval, a strongly compressed single femtosecond pulse of 
multi-millijoule energy and dramatically increased peak intensity will be generated. Input pulses of parameters 
used in these simulations can be generated by the present oscillator-chirped pulse ampUfier Titan:Sapphire laser 
systems, which makes the results of this study suitable for practical applications. 
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Figure 8: Development of the spatiotemporal pulse shape (a "time view") and the respective parameters of the 
pulse taken at position of maximal time compression, z = zsc< versus the gas pressure of the medium. 



All results found here directly follow from the solution of the NLSE. From the other side, the behavior of 
the pulse in regime of SC can be physically explained if we consider the interplay between the most active 
and important processes, in this case - the SF and the dispersion. In the NLSE, they are described by the 
terms ikn2nQ^\E\'^E (in the space domain) and z(/?2/2)9^£'/9r^, respectively. Within such an approach, the 
slower development of the pulse compression in helium than in neon. Figure 1, can be explained with the lower 
nonhnear refractive index 712 ofheliumand, at same peak intensity, smaller nonlinear term zfcn2no^|-Ep.E. This 
makes the process of SF in helium more gradual than in neon. From the other side, the higher value of intensity 
gain, and, as a result - higher value of SC, in helium can be explained by the lower value of GVD of helium. 
While more slowly, the nonhnear term succeeds to SF the pulse in helium at almost same spot diameter. Figure 
1(a). At the same time, the lower GVD in helium prevents the pulse from substantial time spreading. Thus, 
the common action of both terms "squeeze" stronger the pulse energy in the spatiotemporal domain, which 
results in a higher increase in the "vertical" direction, i.e., intensity gain, and, thus - stronger SC, Figure 1. The 
higher intensity gain and pulse compression at low input energy. Figure 2 and Figure 4, also looks naturally if 
it is considered within the interplay between above two terms. The low input energy for pulses of same initial 
pulse duration and same size in transversal direction means low input intensity and field amplitude. The latter 
also means low steepness of the pulse fronts and, thus, relatively low value of the time derivatives. As a final 
result, this diminishes the expanding rate due to the dispersion term i(/32/2)(?^£'/9T^ in which, the second time 
derivative determines dynamically the magnitude of the expanding rate. This explains why at low input energy, 
the magnitude of the SC is higher. The above factors support the opposite trend with the increase of the input 
energy of the pulse. The higher input energy means higher input peak intensity and, consequently, higher value 
of the field time derivatives and higher magnitude of the dispersion term. The latter intensifies the process of the 
dispersion spreading and spUtting, thus preventing the pulse to achieve higher levels of intensity gain and SC 
(in comparison with the lower energy case). The results in Figure 3 and Figure 5 also look transparent within 
the framework of such interplay. At low energy, but equal other conditions, the peak intensity and the field 
ampUtude of the initial pulse is relatively low, which results in low rate of development of the self-focusing by 
the cubic term of the NLSE, ikn2n'^^^\E\'^E. That is why, the position where the maximal SF, and, as a result, 
the maximal intensity gain and time compression occur is far from the entrance point of the medium. Figures 
3,5. Increasing the input pulse energy, and thus, the input peak intensity, intensifies the SF and the pulse needs 
shorter distance, zjq and zsc, to develop the maximal intensity gain and SC, respectively. Finally, the increased 
peak intensity gain. Figure 6, and, as a result, stronger SC (shorter pulse durations). Figure 7, with reduction 
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of gas pressure can be explained by the reduced value of the GVD of the medium, and in this way - reduced 
value of the whole dispersion term at lower pressure. This diminishes the dispersion spreading and promotes 
the intensity gain and the SC. Of course, the reduced pressure reduces also the (total) nonlinear refractive index 
of the medium, and, in this way, the whole nonlinear term. The reduction of the nonlinearity results in slower 
rate of development of the SF and extends distances ziq and zgc to generate pulses of maximal intensity and 
SC, see the data of ziq and z^c in Figures 6 and 7. The numerical simulations show that the net effect of the 
reduced dispersion dominates, and the final effect is increased pulse compression. 

The existence of SC before the splitting has crucial importance for the application of this effect. For example, 
it allows, to certain degree, to overcome, or, at least, to relax the requirement for existing of stable soliton-like 
pulse propagation for some applications. The theoretical and experimental studies show that such pulses do not 
exist within the specified basic set of optical processes and experimental conditions, considered here. For many 
applications, however, it is not necessary to have a stable pulse along the whole propagation distance but only a 
pulse of specified parameters at given space point. Such point will be called space target point. The generation of 
such pulses can be considered as a satisfactory alternative to the soliton generation at conditions not supporting 
soUtons. In this way, the soliton concept (where it is not applicable) can be replaced by a controllable guidance 
concept, which we introduce here. The controllable guidance concept seems to be appUcable in broad range of 
conditions, in which the pulse propagation dynamics can be unambiguously predicted. The absence of ionization 
in our case plays important role in the realization of such a concept because it ensures a smooth and irreversible 
(along the propagation distance of interest) pulse propagation dynamics. The latter consists in a single SC 
event, reaching maximal efficiency at 2:sc> followed by slight broadening of the pulse along a short distance 
(from ^;sc to Ziq), while the peak intensity continuously increases (reaching maximal value at 2; = Ziq), and, 
finally, initiation/developmentof splitting of the pulse at/beyond z = ziq. Figure 1. Although the pulse evolves 
along the whole propagation distance, the SC mechanism together with the smooth pulse propagation dynamics 
facilitates the generation of a pulse of well predictable parameters at a desired space target point. It is based 
on the fact that the magnitude of the time compression and the intensity gain as well as the space position of 
the strongest time "focusing" z^c and the strongest intensity gain ziq, can be controlled by the pulse {e.g., 
energy), Figures 2-5, and the medium (density/pressure), Figures 6, 7, parameters. The pulse shape, Figure 8, 
also changes smoothly and could be a subject to control by the initial parameters. Such a control, however, has 
few shortcomings. At first, it is a passive kind of control and, once launched in the medium, the pulse is no 
more subject to external control. In the second place, such a control is not complete because the parameters of 
the compressed pulse cannot be controlled independently. In any case, however, the investigations toward the 
controllable guidance concept look perspective. 




Pulse Duration [fs] 

Figure 9: Time \E{r = Q,z = zsc, t) |^ intensity profiles of the tsunami pulse at the space position of maximal 
time compression, at P = 17 atm of helium. Sech^ time shape (dashed line) is also shown for comparison. 



A very interesting feature of the pulse propagation can be distinguished from the results in Figures 1 and 
8. As can be seen, the initial pulse has relatively long pulse duration and relatively low peak intensity. The 
parameters of the compressed pulse are most interesting for the practical applications around the points of 
maximal pulse compression, ^sc. and/or the maximal intensity gain, zsc, which play role of target points. Far 
from given space target point, the pulse develops slowly and the peak intensity increases almost linearly with 
the propagation distance. Figure 1(b). Approaching the target point, the dynamical rearrangement of the pulse 
rapidly speeds up and its intensity increases in an exploding manner. Figure 1(b), whereas the pulse duration 
shortens gradually. Figure 1(c). This is also illustrated by Figure 8 if compare the initial pulse and some of 
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the compressed pulse, say, at the lowest pressure - P — 15 atm. Such a behavior of an optical pulse closely 
resembles the behavior of the water tsunami wave, approaching the shore. That is why, such a compressed 
optical pulse whose intensity dramatically surges around the target points will be called optical tsunami pulse. 
The time shape of the intensity profile of the compressed pulse at zsc and P = 17 atm of helium is shown 
in Figure 9. Although the input pulse has Gaussian temporal shape, the intensity time profile of the tsunami 
pulse is very close to sech^ pulse shape, also shown for comparison in Figure 9. The latter represents the 
temporal shape of the fundamental soliton of the (1+1)D NLSE in the case of medium of negative GVD Ii331 . 
Insubstantial deviation of the shape of the tsunami pulse from that one of the fundamental soliton exists only 
around the wings of the pulse. Consequently, in accordance to the controllable guidance concept, a strongly 
compressed femtosecond optical tsunami pulse of huge intensity and soliton-like time intensity profile can be 
generated in given space point and the parameters of that pulse can be controlled within given limits by the input 
pulse and the medium parameters. 

4 Conclusions 

Self-compression of multi-millijoule femtosecond laser pulses has been found as a result of the interplay of the 
diffraction, group velocity dispersion of second order and Kerr nonlinearity of third order in lighter noble gases, 
neon and helium, at conditions of positive group velocity dispersion. These three processes are sufficient to 
trigger the self-compression at a low magnitude of positive group velocity dispersion of the medium. A com- 
parative study of the self-compression in neon and helium has been performed and it confirms the importance of 
the low value of group velocity dispersion for the operation of the self-compression mechanism. The pressur- 
ized helium appears to be the best medium for compression the high energy femtosecond laser pulses within the 
pulse compression concept considered here. The parameters of the compressed pulse can be controlled within 
certain limits by the initial pulse and the medium parameters. A generation of compressed pulse of strongly 
increased intensity and soliton-like temporal profile in the vicinity of given space target point, called optical 
tsunami, is predicted. 
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